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ABSTRACT 

Two massive binary systems of the microquasar type, LS 5039 and LSI +61° 303, 
have been suggested as possible counterparts of EGRET sources. LS 5039 has been also 
recently detected in TeV gamma-rays. Since the massive stars in these binary systems 
are very luminous, it is expected that high energy gamma-rays, if injected relatively 
close to the massive stars, should be strongly absorbed, initiating inverse Compton e 
pair cascades in the anisotropic radiation from stellar surfaces. We investigate influence 
of the propagation effects on the spectral and angular features of the 7-ray spectra 
emerging from these two binary systems by applying the Monte Carlo method. Two 
different hypothesis are considered: isotropic injection of primary gamma-rays with the 
power law spectrum due to e.g. interaction of hadrons with the matter of the wind, 
and the isotropic injection of electrons, e.g. accelerated in the jet, which comptonizc 
the radiation from the massive star. It is concluded that propagation effects of 7-rays 
can be responsible for the spectral features observed from LS 5039 (e.g. the shape 
of the spectrum in the GeV and TeV energy ranges and their relative luminosities). 
The cascade processes occurring inside these binary systems significantly reduce the 
gamma-ray opacity obtained in other works by simple calculations of the escape of 
gamma-rays from the radiation fields of the massive stars. Both systems provide very 
similar conditions for the TeV 7-ray production at the periastron passage. Any TeV 
7-ray flux at the apastron passage in LSI +61° 303 will be relatively stronger with 
respect to its GeV flux than in LS 5039. If 7-rays are produced inside these binaries 
not far from the massive stars, i.e. within a few stellar radii, then clear anticorrelation 
between the GeV and TeV emission should be observed, provided that primary 7-rays 
at GeV and TeV energies are produced in the same process by the same population of 
relativistic particles. These 7-ray propagation features can be tested in the near future 
by the multi-wavelength campaigns engaging the AGILE and GLAST telescopes (>30 
MeV) and the Cherenkov telescopes (>100 GeV, e.g. MAGIC, HESS, VERITAS and 
CANGAROO). 
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1 INTRODUCTION 

Massive binary systems provide very promising conditions 
for acceleration of particles to relativistic energies and also 
well defined conditions for their possible interaction. There- 
fore, they have been considered for a long time as possible 
sources of high energy 7-rays and neutrinos in which accel- 
eration processes can be tested. In fact, some GeV 7-ray 
sources observed by EGRET detector (> 100 MeV) have 
been proposed to be related to well known massive bina- 
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ries in which non-thermal processes were evident at lower 
energies, e.g. LSI +61° 303 (2EG J0241+6119; Thompson 
et al. 1995), Cen X-3 (Vestrand, Sreekumar & Mori 1997), 
Cyg X-3 (2EG J2033-4112; Mori et al. 1997), and LS 5039 
(3EG J1824-1514; Paredes et al. 2000). Early searches of the 
TeV 7-ray signals from massive binaries have not been very 
convincing, with the exception of Cen X-3 binary system 
(containing slowly rotating neutron star) which has been re- 
ported as TeV 7-ray source by the Durham group (Chadwick 
et al. 1998, 1999, Atoyan et al. 2002). The turning-point 
came recently with the observations of TeV 7-ray signals 
from two massive binaries: PSR B1259-63/SS 2883, contain- 
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ing a radio pulsar with the period of 47.8 ms (Aharonian et 
al. 2005a), and LS 5039, so called microquasar possibly con- 
taining a solar mass black hole (Aharonian et al. 2005b). 

The 7-ray emission from massive binaries is usually in- 
terpreted in terms of the inverse Compton scattering (ICS) 
model in which thermal radiation coming from the stellar 
surface is scattered by electrons accelerated in the pulsar 
wind shock (e.g. Maraschi & Treves 1981) or by electrons 
moving highly anisotropically in the form of beams or jets 
(e.g. Bednarek et al. 1990). The binary systems which are 
considered as an example in this paper, LS 5039 and LSI 
+61° 303, have been recently discussed in terms of the mi- 
croquasar IC model by Bosch-Ramon & Paredes (2004a,b) 
and Paredes, Bosch- Ramon & Romero (2005). Primary 7- 
rays might be also produced in the microquasar scenario of 
the binary system as a result of the interaction of hadrons, 
accelerated in the jet, with the matter of the stellar wind 
(Romero, Christiansen & Orellana 2005). For review of the 
7-ray production in microquasars see e.g. Romero (2004) or 
Paredes (2005). 

The problem of importance of the propagation of TeV 
7-rays inside compact massive binary systems appeared at 
the early 80's after first reports on possible observation of 
such 7-rays (see e.g. Weekes 1988). The optical depths for 
7-rays in the radiation field of the accretion disk around the 
compact object inside the binary system have been calcu- 
lated by e.g. Carraminana 1992 and Bednarek (1993) and in 
the radiation field of the massive star by e.g. Protheroe & 
Stanev 1987 and Moskalenko, Karakula & Tkaczyk 1993). In 
the case of very compact binaries, the TeV 7-rays injected 
close to the surface of the massive star initiate the Inverse 
Compton e ± pair cascade. The conditions for which the cas- 
cade processes should become important have been consid- 
ered by Bednarek (1997, Sect. 2) and Dubus (2005, Sect 6). 
In general, the product of the massive star luminosity and 
the square of its surface temperature should be much greater 
than ~ 10 45 erg s _1 K 2 (see Eg. 1 and Fig. 3 in Bednarek 
1997). The 7-ray spectra which emerge toward the observer 
from such compact binary systems strongly depend on the 
phase of the injection place of the primary relativistic par- 
ticles in respect to the observer. At some phases significant 
TeV 7-ray fluxes are expected but at others only photons 
with energies extending to a few tens GeV are able to es- 
cape (see detailed calculations of such anisotropic cascades 
and their application to specific sources, e.g. Cen X-3 and 
Cyg X-3, in Bednarek 1997 (B97), Bednarek 2000 (BOO), 
and Sierpowska & Bednarek 2005 (SB05)). Recently, optical 
depths of TeV 7-rays have been calculated for other TeV 
7-ray sources, LS 5039, PSR 1259-63 and LSI +61° 303, 
without taking into considerations the effects of the IC e ± 
pair cascading (Bottcher & Dermer 2005 and Dubus 2005). 

Note that only calculations by Bednarek and collabora- 
tions (B97, BOO, SB05) and Dubus (2005) take into account 
dimensions of the massive star. Therefore, they can be ap- 
plied to very compact binaries in which injection distance 
of the TeV 7-rays from the center of massive stars is be 
comparable to their radii, e.g. at the periastron passage of 
the compact objects in LS 5039 and LSI +61° 303, or in the 
case of propagation of 7-rays injected at larger distances but 
passing close to the surface of the massive stars. The com- 
parison of the exact calculations with the approximate ones, 




Figure 1. Schematic picture of a compact binary system com- 
posed of a massive star and a compact object (a black hole or 
an neutron star) on an orbit around the massive OB star. The 
observer is located at the inclination angle 9 and the azimuthal 
angle u>. The matter accreting onto a compact object from the 
massive star creates an accretion disk. Particles (electrons or pro- 
tons) are accelerated inside the jet launched from the inner part 
of an accretion disk. Primary electrons and/or 7-rays, injected 
at the distance z from the base of the jet, initiate an anisotropic 
inverse Compton pair cascade in the radiation field of the 
massive star. A part of the primary 7-rays and secondary cascade 
7-rays escape from the binary system toward the observer. 



which neglect dimensions of the massive stars, can be found 
in Dubus (2005). 

In this paper we concentrate on details of the propa- 
gation of high energy 7-rays injected close to the surface of 
the massive star taking into account the effects of cascades 
initiated through the Inverse Compton and e^ pair produc- 
tion processes. We calculate the 7-ray spectra emerging to 
the observer from such anisotropic cascades for two compact 
massive binaries LS 5039 (already reported in GeV and TeV 
7-rays) and LSI +61° 303 (reported only at GeV energies). 
Specific cases with the injection of primary 7-rays (e.g. by 
relativistic hadrons) or electrons (accelerated in the jet or 
the shock) with simple power law spectra and spectral in- 
dexes equal to 2 are considered in order to better understand 
the basic features of such anisotropic cascade processes. 



2 THE MASSIVE BINARY SYSTEMS 

Both binary systems considered in this paper belong to the 
class of the non-thermal radio high mass X-ray binaries 
showing evidences of collimated relativistic outflows. They 
are called microquasars due to their supposed similarities to 
quasars which show very narrow jets moving with relativistic 
speeds. In binary systems jets are launched from the inner 
parts of accretion disks around compact objects (a neutron 
stars or a solar mass black holes). 

Here we consider a simple scenario in which jets are 
launched along the disk axis. The surface of the disk is in 
the plane of the binary system, i.e. the jet direction is per- 
pendicular to the plane of the binary system. The schematic 
picture of such binary system is drawn in Fig.0 Relativistic 
particles, injected in the jet at the distance z from its base, 
produce 7-rays. If the injection site of the 7-rays is relatively 
close to the massive star, 7-rays interact with the stellar ra- 
diation initiating the IC e pair cascade. The efficiency of 
such cascades depend strongly on the parameters of the bi- 
nary systems LS 5039 and LSI +61° 303. These two binaries 
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differ in basic parameters, e.g their orbital periods are equal 
to 3.9 days and 26.5 days. However both of them provide 
conditions in which cascading effects have to be taken into 
account when considering 7-ray production processes. 



2.1 The binary system LS 5039 

This binary system shows relativistic radio jets on miliarc- 
second scales, with the speed of v ~ 0.3c (Paredes et al. 
2000). It has been also suggested to be a counterpart of 
the EGRET source 3EG J1824-1514 localized at ~ 0.5° 
(Paredes et al. 2005). This source has relatively flat spec- 
trum above 100 MeV, spectral index < 2, and the 7-ray 
luminosity ~ 4 x 10 35 erg s _1 . Moreover, the position of 
LS 5039 is consistent at the 3cr level with recently detected 
TeV source HESS J 1826- 148 (Aharonian et al. 2005). The 
spectrum above 250 GeV is also flat with the photon in- 
dex 2.12 ± 0.15, although the luminosity is only ~ 10 33 erg 
s~ , about two orders of magnitude less than at GeV ener- 
gies. Recent analysis of the TeV 7-ray light curve by Casares 
et al. (2005b), using new orbital parameters, show possible 
flux variations of a factor ~ 3 with the maximum around 
the phase ~ 0.9. 

The basic parameters of the binary system LS 5039 has 
been recently reported by Casares et al. (2005b): the semi- 
major axis, a = 3.4r*, ellipticity e = 0.35 ± 0.04, the az- 
imuthal angle of the observer in respect to the periastron 
passage, lj = 225°, radius of the massive star, = 9.3lQg 
Rq, and its surface temperature, T B = 3.9 x 10 4 K. For 
these parameters the distance of the compact object from 
the massive star chances in the range: r p = 2.2r* at the 
periastron up to r a = 4.5r* at the apastron. The estimated 
inclination angle of the binary system toward the observer 
depends on the mass of the compact object. It is estimated 
on 9 = 24.9 ± 2.8° for the case of the black hole with the 
mass 3.7M0 and on ~ 60° for the neutron star (Casares et 
al. 2005b). We present the results of numerical calculations 
for both inclination angles in the case of LS 5039. 

2.2 The binary system LSI +61° 303 

LSI +61° 303 has been also observed as a non-spherical ra- 
dio source, which structure was interpreted as due to rela- 
tivistic radio jets, the speed of ~ 0.6c, with some hints of 
its precession (Massi et al. 2004). It has been pointed out 
(Gregory & Taylor 1978) that this source is connected with 
the COS B 7-ray source CG135+01 (Hermsen et al. 1977). 
CG135+01 has been also detected by EGRET (MOO MeV); 
the source 2EG J0241+6119 has a hard spectrum with pho- 
ton spectral index 2.05 ± 0.06 (Kniffen et al. 1997), and by 
COMPTEL in the energy range ~0.75 to 30 MeV, spectral 
index 1. 951^3 (van Dijk et al. 1996). The analysis of dif- 
ferent EGRET observations shows evidences of variability 
(Tavani et al. 1998, confirmed by Wallace et al. 2000) with 
probable modulation with the orbital period of LSI +61° 
303 and the maximum emission near the periastron passage 
(Massi 2004). The Whipple group (Hall et al. 2003) puts an 
upper limit on the TeV flux from this source, ~ 10" 11 cm" 2 
s" 1 above 500 GeV (< 1.3 x 10 32 erg s" 1 ), which is clearly 
below an extrapolation of the EGRET spectrum. 

The basic parameters of the binary system LSI +61° 



303 has been recently reported by Casares et al. (2005a): 
the semimajor axis a — 5.3r*, ellipticity e = 0.72, the in- 
clination of the binary system toward the observer is not 
well constrained by the observations (Casares et al. 2005): 
25° < i < 60° for a neutron star and 9 < 25° for a black 
hole. We apply the value of 9 = 30°. The azimuthal angle of 
the observer in respect to the periastron passage is uj = 70°, 
radius of the massive star r* = 13.4 Rq, and its surface 
temperature T B — 2.8 x 10 4 K. For these parameters the dis- 
tance of the compact object from the massive star changes 
in the range: r p = 1.5r* at the periastron up to r a = 9.15r* 
at the apastron. We apply these parameters in our further 
calculations. 



3 OPTICAL DEPTHS FOR 7-RAYS 

The optical depths for 7-rays (for the process 7 + 7 — > e + e~) 
in the radiation field of the massive star, determined by the 
radius of the star and its surface temperature, are calculated 
in the general case, i.e. for an arbitrary place of injection of 
7-ray photons with arbitrary energies and angles of propa- 
gation (see also B97, BOO, SB05). This approach allows us 
to calculate the optical depths even for the primary 7-rays 
injected at the surface of the massive star. Therefore, it can 
be applied for studies of the cascade processes initiated by 
7-ray photons since secondary 7-rays can appear, in princi- 
ple, everywhere inside the binary system. In fact, the optical 
depths for the 7-rays propagating in the thermal radiation 
of these massive stars can be obtained by simple re-scaling 
of the earlier calculations for the massive star in Cen X-3 
(shown e.g. in Bednarek 2000, Fig. 2), since they are propor- 
tional to the fourth power of the surface temperature and 
the square of the radius of the massive star. Also a shift 
in energy proportional to the surface temperature is neces- 
sary. Note, that in the case of cascade processes discussed 
in this paper, primary particles propagating at specific di- 
rection can contribute to the final 7-ray spectra escaping 
at other directions. Therefore, for easier analysis of the ob- 
tained results, we show here the optical depths for 7-rays 
injected at the distance of the periastron and the apastron 
passages of the compact object around the massive stars 
in LS 5039 (Fig. Hi,b) and LSI +61° 303 (Fig. EJj,d), as a 
function of the photon energy and their arbitrary injection 
angles, a, measured from direction defined by the centers of 
the companion stars, toward the massive star. 

As expected, the optical depths strongly depend on the 
injection parameters (photon energy, angles of propagation) 
and the parameters of the massive star. The optical depths 
reaches the maximum corresponding to the peak in the black 
body spectrum of soft photons. The maximum optical depth 
shifts to larger photon energies with decreasing angle a (pro- 
vided that q is smaller than the angle, (5, intercepted by 
the massive star observed from the distance of the injec- 
tion place). The optical depths are lower for a > n — [3, 
since in this case 7-rays propagate only to the surface of 
the massive star (see full thick curves in Figs. |2K and c). 
7-rays with energies ~ 1 TeV, injected at the periastron 
distance in LS 5039, 2.2r*, are absorbed (r > 1) for most 
of the propagation directions, except a small cone with the 
angular extend of ~ 40°. At the apastron distance, 4.5r*, 
the escape cone increases to ~ 60°. Only within such small 
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Figure 2. The optical depths for 7-rays (as a function of their energy) on e pair production in collisions with stellar photons. 7-rays 
are injected at two distances from the massive star in LS 5039 (2.2r* (a) and 4.5r* (b)) and in LSI +61° 303 (1.5r* (c) and 9.15r* (d)), 
corresponding to the periastron and the apastron passages, respectively. Specific curves show the optical depths for the injection angles 
of 7-rays measured from the direction defined by the centers of the stars: a = 150° (full curve, direction toward the massive star), 120° 
(dashed), 90° (dotted), 60° (dot-dashed), and 30° (dot-dot-dot-dashed). The optical depths are calculated up to the infinity accept two 
cases (marked by full curves in (a) and (c)) which are calculated to the surface of the massive star. 



cones, 7-rays with ~ 1 TeV energy have high chance to es- 
cape without absorption. The optical depth for 7-rays in 
LSI +61° 303 are typically lower due to significantly lower 
surface temperature of the massive star. At the periastron 
distance, 1.5r*, the escape cone for ~ 1 TeV 7-rays is also 
~ 40°, very similar to the periastron distance in LS 5039, 
due to the closer location of the injection place. However, 
at the apastron distance, 9.157%, most of ~ 1 TeV 7-rays 
escape from the binary system without absorption, except 
7-rays moving toward the massive star within the cone with 
the angle ~ 80°. 

In LS 5039, the angles between the observer (at the in- 
clination angle 25°) and the direction defined by the stars 
are ~ 110° for the periastron passage and ~ 70° for the apas- 
tron passage. For such geometry, the optical depths towards 
the observer are larger than unity for 7-rays with energies in 
the range ~ 0.03 — 20 TeV, for the periastron passage, and 
~ 0.2 - 2 TeV, for the apastron (see Fig.|5Ji,b). These cal- 
culations of the optical depths are generally consistent with 
that obtained by Dubus (2005). 

In LSI +61° 303, the observer is located at the angle 
of ~ 80° at the periastron passage and at ~ 100° at the 
apastron passage. In this case the optical depths are larger 
than unity for 7-rays with energies in the range ~ 0.1 — 10 
TeV for the periastron passage, but always lower than unity 
for the apastron, provided that they are injected toward the 
observer (see Fig. |5J;,d). 



4 THE CASCADE 7-RAY SPECTRA 

Since the optical depths for high energy 7-rays in the radi- 
ation fields of the massive stars in binary systems LS 5039 
and LSI +61° 303 are large, the total 7-ray spectra, emerg- 
ing from the binary systems toward the observer, are deter- 
mined by the primary 7-ray spectra (if such are produced 
e.g. by hadrons accelerated in the jet) and by the 7-ray 
spectra produced in the IC e pair cascades occurring in 
the radiation of the massive star. However, since the injec- 
tion place of primary 7-rays or primary electrons is different 
than the center of the massive star (the source of isotropic 
soft radiation), the IC e± pair cascade occurs in fact in the 
anisotropic radiation field. IC e pair cascades of such type 
have been already considered by us under two extreme as- 
sumptions. In the first approximation, local isotropization of 
the secondary cascade e ± pairs by the random component 
of the magnetic field inside the stellar wind and the jet is 
assumed (see B97, BOO). In the second approach, we follow 
the paths of secondary cascade e pairs in the magnetic field 
which structure is described by a specific model (see SB05) . 
In the present work we apply the first approach. Note, that 
the IC e pair cascade considered in this work does not take 
into account the synchrotron losses of e pairs and assumes 
that secondary e pairs radiate locally secondary 7-rays, 
i.e. at their production site. The conditions for such type 
of cascade are determined in Bednarek (B97, see Sect. 2). 
For completeness and due to some different details of the 
scenario considered in this paper, we discuss below some 
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important conditions for IC e pair cascades occurring in 
the radiation field of the massive star. 

Leptons injected into the radiation of the massive star 
(primary electrons and secondary e pairs from the cascade) 
are immersed in the stellar wind and/or in the plasma out- 
flow along the jet. Therefore, efficiency of the ICS process 
is determined by the relative importance of the IC cooling 
time scale in respect to the characteristic escape time scale, 
determined by the advection time scale of e pairs with the 
stellar and/or jet outflows. The IC cooling time of leptons, 
the Thomson regime, in the radiation of the massive star 
can be estimated from, 

ric = m c7c /Pf n c « 4 x 10V/( 7c T 4 4 ) s, (1) 

where m c and 7 e are the rest mass and the Lorentz factor of 
leptons, P^ = (4/3)c<tt7c is the energy loss rate on ICS of 
leptons, c is the velocity of light, <tt is the Thomson cross 
section, [7 ra< j = ctsbTb /r 2 , <tsb is the Stefan-Boltzmann con- 
stant, T s = 10 4 T4 K and R — rr* are the surface tempera- 
ture of the massive star and the distance from the center of 
the massive star (in units of its radius). The characteristic 
escape time scale of leptons from their creation (accelera- 
tion) place, identified with their advection time scale with 
the jet or the wind plasma flow, is 

r csc = R/V « 33r*,i 2 r/« s, (2) 

where v = V/c is the velocity of the jet (or the stellar wind). 
It is assumed that typical radii of the massive stars are of 
the order of r* = 10 12 rv,i2 cm. By comparing Eq. 1 and 
Eq. 2, we estimate the minimum Lorentz factor of leptons 
above which they cool locally, 

7c > 7 min « 1.2 x 10 4 r«/(T 4 V,i 2 ). (3) 

For example, in the case of the massive star in LS 5039, 
T 4 = 3.9, r*,i2 = 1, and v = 0.5, leptons with the Lorentz 
factors above, 7m i n ~ 25r, cool locally before escaping from 
the binary system. 

Leptons with Lorentz factors, 7 m in, produce 7 -ray pho- 
tons with energies, 

E- t « (4/3)e 7 L n « 10 3 rV/m 7 < 12 ) MeV, (4) 

where typical energies of photons coming from the massive 
stars are e = 3k B T B « 2.6 x 10~ 6 T 4 MeV. We conclude that 
spectra of 7-ray photons produced in the cascade (assuming 
local cooling of leptons) are correct above energies given by 
Eq. 0| For the parameters of the massive stars in LS 5039 
and LSI +61° 303, escaping 7-ray spectra are correct above 
~ 10 3 MeV, provided that production of 7-rays occur inside 
the jet and within r ~ 30r* from the center of the massive 
star. For the stellar wind region, the limits on the low energy 
cut-offs in the cascade 7-ray spectra are about two orders of 
magnitude lower due to much lower stellar wind velocities 
in respect to the plasma velocity in the jet. 

However, for the Lorentz factors above ~ 10 scattering 
of soft photons occurs in the Klein-Nishina (KN) regime. We 
approximate the IC cooling time of leptons in the KN regime 
by, 

ric N = m c7c /P? ? c (7KN/T) « 8 x 10-V 7o /r 4 2 s, (5) 

where 7 kn/t ~ m c c 2 /e ~ 2 x 10 5 /T4. By comparing the 
escape time scale, r esc , with the IC cooling time scale in the 



KN regime r IC , we obtain the upper limit on the Lorentz 
factor of leptons which are able to cool locally 

7c < 7max ^ 4 x 10V,i 2 T 4 2 /(™). (6) 

Therefore, we conclude that leptons with energies ~ 10 TeV 
should be able to cool locally within a few stellar radii from 
the surface of the massive stars in LS 5039 and LSI +61° 
303. 

In the calculations shown below we assume that syn- 
chrotron energy losses of leptons can be neglected in respect 
to the IC losses. This is correct provided that the surface 
magnetic fields of the massive stars fulfill the condition, 

B s < 40T 4 2 G, (7) 

obtained from the comparison of the synchrotron and IC 
energy losses of leptons in the Thomson regime. The sim- 
ple limit given by Eq. Q is valid under assumption on the 
B oc r -2 magnetic field dependence on the distance from the 
stellar surface. It neglects the inner dipole part of the mag- 
netic field in which B oc r~ 3 . Therefore, real upper bound on 
the surface magnetic field is a factor of 2-3 larger than given 
by Eq.0 For the massive stars in LS 5039 and LSI +61° 303, 
these upper bounds (Eq. are ~ 620 G and ~ 320 G. In 
the Klein-Nishina regime, i.e. for 7 e S> 10 5 , the limit on the 
magnetic field is more restrictive due to the decrease of the 
cross section. It has been discussed in details by Bednarek 
(B97, see Fig. 4 in that paper). 

Under conditions specified above (neglected syn- 
chrotron losses of leptons and their local isotropization) , we 
study the features of the 7-ray spectra emerging from the 
binary system, due to the propagation effects in the radi- 
ation of their massive companions, assuming that primary 
particles, i.e. electrons or 7-rays, are isotropically injected 
somewhere along the jet. We consider the primary particles 
with the power law spectra and spectral index equal to 2 
due to the equally distributed power per decade. The high 
energy cut-off in these primary spectra is assumed at 10 
TeV (to be consistent with the observations of 7-ray spec- 
trum from LS 5039 up to ~ 4 TeV, Aharonian et al. (2005)). 
These initial spectra of particles are normalized to unity. The 
spectral index equal to 2 is also motivated by a relatively flat 
GeV 7 -ray spectra observed from these two sources. The jets 
in microquasars move with much lower velocities (Doppler 
factor D ~ 1) than these ones observed in active galactic 
nuclei. Therefore, the effects of relativistic beaming can be 
neglected in the first approximation. 



4.1 LS 5039 

Let's at first consider the case of isotropic injection of pri- 
mary electrons at the base of the jet, i.e. z m 0. In fact, 
7-rays have to be injected at some distance from the accre- 
tion disk at which the disk radiation can be neglected in 
respect to the stellar radiation. The distance from the base 
of the jet at which the above condition can be fulfilled is 
estimated by comparing the energy density of the disk radi- 
ation with the energy density of the stellar radiation. These 
radiation fields are defined by the surface temperature of 
the massive star and the maximum temperature on the sur- 
face of the disk (at the disk inner radius, Tl n ). We approxi- 
mate the accretion disk radiation by the model of Shakura 
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Figure 3. Differential 7-ray spectra (multiplied by the square of photon energy) escaping from the binary system at a specific range of 
the cosine angles a, measured in respect to the direction defined by the injection place and the center of the massive star. The range of 
A cos a (with the width 0.1) are centered on 0.95 (full curve), 0.55 (dashed), 0.15 (dot-dashed), -0.25 (dotted), and -0.65 (dot-dot-dot- 
dashed). 7-rays are produced in the cascade initiated by primary 7-rays with the power law spectrum and spectral index 2 which are 
injected isotropically close to the base of the jet (z ss <C r*, but sufficiently far away from the accretion disk), and at two distances 
from the massive star 2.2r* (upper figures, corresponding to the periastron passage of the compact object) and 4.5r* (bottom figures, 
corresponding to the apastron passage), (a) Spectra of primary 7-rays escaping without interaction, (b) Spectra of secondary 7-rays 
produced in the IC cascade, (c) Total spectra of 7-rays escaping from the binary system (the sum of spectra shown in (a) and (b)). 
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Figure 4. As in Fig. |3]but for the injection distance along the jet z = 5r*. The distance of the injection place from the massive star is 
~ 5.5r* at the periastron (upper figures) and ~ 6.7r* at the apastron (bottom figures). 
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Figure 5. Total 7-ray spectra from the cascade initiated by isotropically injected primary leptons, with the power law spectrum and 
spectral index 2, at the range of cosine of the observation angles cos a, with the width A cos a = 0.1. The spectra are shown for cos a 
centered on 0.95 (full curve), 0.55 (dashed), 0.15 (dot-dashed), -0.25 (dotted), and -0.65 (dot-dot-dot-dashed), for two distances of the 
injection place from the massive star 2.2r* (figures (a) and (c)) and 4.5r* (figures (b) and (d)), corresponding to the periastron and the 
apastron of the injection place, respectively. The spectra are shown for the injection place at the base of the jet z f» <C ((a) and (b)) 
and at the distance of z=5r* from the base of the jet ((c) and (d)). 



& Sunyaev (1973) in which the disk radiation from its sur- 
face can be approximated by the black body radiation with 
some temperature gradient. The following condition has to 
be approximately fulfilled, 

T s 4 /r 2 « H 4 n (r in /z) 2 , (8) 

where n n is the disk inner radius. From this condition, we 
estimate that above, 

z ~ r in r(T in /T s ) 2 , (9) 

the massive star radiation dominates over the accretion disk 
radiation. For typical parameters of the considered binary 
systems, T s = 3 x 10 4 K, r — 10, the accretion disk in- 
ner radius n n = 10 7 cm, and T ln — 10 6 K (limited by the 
condition that the disk thermal luminosity has to be lower 
than the observed X-ray luminosity from these binary sys- 
tems, < 10 35 erg s _1 ), we estimate that above z ~ 10 11 cm 
(~ O.lr*) from the base of the jet, stellar radiation domi- 
nates over accretion disk radiation. This condition is valid 
in the Thomson regime. However for the IC scattering pro- 
cess occurring in the Klein-Nishina regime, stellar radiation 
starts to dominate even at lower distances from the inner 
part of the disk, due to significantly larger temperature in 
the inner disk than on the stellar surface. These same argu- 
ments concern also possible disk corona which energy den- 
sity can be comparable to the energy density of radiation 
from the surface of the disk but the characteristic temper- 
ature (average photon energies) are larger. The conditions 
for absorption of 7-ray photons in the radiation of the mas- 
sive star are more favorable in respect to the radiation of 



the accretion disk due to larger angular extend of the stellar 
disk in respect to the inner part of the accretion disk. In 
conclusion, we consider the injection of the primary 7-rays 
and electrons at the base of the jet assuming that it occurs 
at the distance, z > O.lr* w 0, which is very close to the 
base of the jet in respect to dimensions of the massive star. 

Although primary 7-rays and soft photons from the 
massive star are injected isotropically, their injection places 
are located at different parts of the binary system (the jet or 
the compact object for primary 7-rays and the massive star 
for soft photons). Therefore, in fact primary 7-rays develop 
an IC e* pair cascade in the non-isotropic radiation field. 
The spectra, which emerge from the binary system, depend 
on the location of the observer in respect to direction de- 
fined by the location of the injection place and the center 
of the massive star (the angle a measured from the outward 
direction defined by the injection place and the center of the 
massive star). We calculate the angle dependent spectra of 7- 
rays escaping from the binary system, by applying the Monte 
Carlo method (details are described in Bednarek BOO). Such 
method allows us to include the re-distribution of directions 
of secondary cascade 7-rays in respect to directions of their 
parent primary 7-rays, which is due to the isotropization of 
secondary cascade e pairs and their preferable head on in- 
teractions with soft photons arriving from specific direction 
on the sky (i.e. inside the stellar disk limb). The parameters 
of 7-ray photons produced in the cascade (their energies and 
the escape angles a) are sorted at a specific range of the co- 
sine angles, cos a, with the width of A cos a = 0.1. 
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We show the spectra of primary 7-rays which escape 
from the radiation field of the massive star without absorp- 
tion (Fig. Et-) t the spectra of 7-rays produced as a secondary 
cascade products (Fig. [SJa) , and the sum of these two, i.e. the 
total 7-rays spectra (Fig.0), for primary 7-rays injected at 
the distance 2 » along the jet, and for two locations of the 
compact object on its orbit, at the periastron distance from 
the massive star (2.2r*, upper figures in Fig. 01 and at the 
apastron distance (4.5r*, bottom figures). Simple absorption 
effects of primary 7-rays, for both distances of the injection 
place from the massive star, are very strong (see Fig. 0i). 
Primary 7-rays with energies in the range between a few tens 
of GeV up to a few TeV are completely absorbed (the exact 
range depends strongly on the observation angles) , provided 
that the observation angles lay within the hemisphere con- 
taining the massive star. However, these strong deficits of 
7-rays (between 0.1 — 1 TeV) are partially fulfilled by the 
secondary cascade 7-rays (see Fig. 0>). The secondary 7- 
rays contribute also to the total 7-ray spectrum below ~ 10 
GeV escaping in the outward directions, i.e. for a > 90° 
(Fig. 0d) . Total 7-ray spectra (Fig. 0) strongly depend on 
the observation angle, they show strong deficit above ~ 100 
GeV (a dip up to two orders of magnitudes) and the excess 
(up to a factor of two) in respect to the shape of the primary 
power law spectrum for large angles a. However, they be- 
come more similar to the injected 7-ray spectrum for small 
a. As expected, 7-ray spectra escaping from the binary sys- 
tem for the injection places located at larger distances from 
the massive star (at the apastron passage), are less mod- 
ified by the cascading processes than these ones produced 
at the periastron passage. The secondary 7-ray spectra pro- 
duced by primary 7-rays injected at the apastron passage 
show decline above a few TeV (Fig. 3b). This is due to the 
high energy cut-off in the spectrum of primary particles at 
10 TeV and due to the average larger interaction angles be- 
tween primary 7-rays and soft photons from the massive star 
at the apastron passage. 

In Fig. we also show the 7-ray spectra produced in 
the cascade for the periastron and the apastron distances 
assuming that the injection place is located at the distance 
of z = 5r* from the base of the jet in order to have impres- 
sion how significant are the 7-ray spectra on the production 
site in the jet. General features of these 7-ray spectra are 
quite similar to the case of injection at the base of the jet 
(z « 0). Therefore, if the injection of primary 7-rays occurs 
with similar acceleration efficiency and spectrum along the 
jet within a few stellar radii from the base of the jet, then 
the angular distribution of 7-rays (and their spectra) formed 
in the cascade process show quite similar features in relation 
to the direction defined by the injection place and the mas- 
sive star. Note however, that the 7-ray spectra toward the 
observer located at fixed direction in respect to the plane 
of the binary may look very different for these two injection 
places as we discuss in Sect. 5. 

As a second scenario we consider the injection of pri- 
mary electrons in the specific region of the jet. As above, we 
consider two locations for the injection place along the jet: 
the base of the jet (Fig. 0i,b) and the distance z = 5r* from 
the base of the jet (Fig. 0,d). The case of isotropic injec- 
tion of primary electrons with the power law spectrum differs 
from previously considered case of isotropic injection of pri- 
mary 7-rays since electrons have tendency of more frequent 



production of the first generation of cascade 7-rays in the di- 
rection toward the massive star (the effect of non-isotropic 
radiation and kinematics of IC process). Therefore, 7-ray 
spectra which escape from the binary system in this sec- 
ond scenario resemble the secondary cascade 7-ray spectra 
produced in the case of injection of primary 7-rays (compare 
Figs. 03 and0> with Fig. 01. Escaping 7-ray spectra strongly 
depend on the viewing angle a not only in TeV range (> 100 
GeV) but also in the GeV range (see Fig. . Moreover, the 
7-ray fluxes expected in these two energy ranges are strongly 
anticorrelated (large TeV flux is accompanied by low GeV 
flux and the opposite). This is due to the fact that at direc- 
tions of the large optical depths the energy converted from 
primary electrons in the IC process to the TeV 7-rays is 
re-distributed to the GeV energies. In directions where the 
optical depths are low, the energy of TeV 7-rays is not so 
efficiently degraded in the cascade process to lower energies. 
As a result, relatively large TeV fluxes are accompanied by 
small GeV fluxes. 

4.2 LSI +61° 303 

The binary system LSI +61° 303 differs in some aspects from 
LS 5039. Although both binaries contain massive stars with 
similar radii, the surface temperature of the star in LSI +61° 
303 is a factor of ~ 0.7 lower which results in the radiation 
energy density by a factor of ~ 4 lower. Also the distance of 
the compact object from the massive star in LSI +61° 303 
change in larger range, from 1.5r* up to 9.15r*. Therefore, 
primary 7-rays injected in the jet propagate through the 
radiation field which vary with larger amplitude with full 
orbital period than in LS 5039. 7-ray spectra escaping to- 
ward the observer for the case of injection of primary 7-rays 
and electrons change more significantly with the observation 
angle a. As an example, we show the 7-ray spectra for the 
case of injection of primary 7-rays at the base of the jet, 
i.e. for z — 0, (Fig. 0, and at the distance z = 3r* from 
the base of the jet (Fig. 01. 7-ray spectra for the periastron 
passage of the injection place are quite similar to these ones 
expected from LS 5039 (compare the results in the upper 
Figs.0 and0 with the upper Figs.0 and 0;). The effects 
of the lower surface temperature of the massive star in LSI 
+61° 303 and the distance at the periastron passage almost 
compensate. However, at the apastron passage and at the 
base of the jet, the absorption dips (above ~ 100 GeV) are 
much lower in the case of LSI +61° 303 (compare the bottom 
Figs. and0 with the bottom Figs. and0). Even at 
the apastron passage of the compact object in LSI +61° 303 
and at the distance of injection place from the base of the 
jet z = 3r*, 7-ray spectrum above ~ 100 GeV should still 
vary by a factor of ~ 5 with the angle a (see Fig. 0) . Note 
that secondary cascade 7-ray spectra produced by primary 
7-rays, which are injected at large distance from the massive 
star (at the apastron passage and at z = 3r*), show strong 
cut-offs after ~ 2 TeV. This is the result of a relatively weak 
absorption of primary 7-rays with energies above a few TeV, 
combined with the kinematics of the IC scattering process 
for the case of strongly anisotropic radiation field of soft 
photons from the massive star. 

7-ray spectra emerging from the binary system LSI 
+61° 303 in the case of injection of primary electrons at the 
periastron distance are very similar to the spectra expected 




log (E / MeV) log (E / MeV) log (E / MeV) 

Figure 6. Differential 7-ray spectra escaping from the binary system at a specific range of the cosine angle a (as in Fig.0 but for the 
binary system LSI +61° 303. The injection place of primary 7-rays is at the distance of the periastron passage (R = 1.5r*, upper figures) 
and the apastron passage (9.15?%, bottom figures) close to the base of the jet (2 fs < r*). The primary 7-rays spectra which escape 
without absorption are in (a), secondary cascade 7-ray spectra in (b), and the sum of (a) and (b), total escaping 7-ray spectra are in (c). 




Figure 7. The 7-ray spectra produced by primary 7-rays (as in Fig. |5J but for the injection place of primary 7-rays and the distance 
z = 3r* from the base of the jet. The distance of the injection place from the massive star is ~ 3.35r s t a r at the periastron (upper figures) 
and ~ 9.6rv at the apastron (bottom figures). 
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Figure 8. As in Fig, 151 but for LSI +61° 303. 7-ray spectra are calculated for the injection place of primary electrons at the periastron 
distance (i.e. at 1.5r* - figures (a) and (c)) and at the apastron distance (i.e. at 9, 15r* - figures (b) and (d)), and at the distance from 
the base of the jet z = (figures (a) and (b)), and z = 3-R* (figures (c) and (d)). 



from LS 5039 (Fig.EJi,c with Fig.0i,c). The absorption fea- 
tures are only slightly stronger. 7-ray spectra, produced by 
electrons injected farther from the base of the jet, e.g. at the 
distance z = 3r*, are still strongly influenced by the cascade 
effects showing strong dependence on the observation angle 
in the GeV and TeV energy ranges (see Fig. [HJi) . However, 
at the apastron distance, the 7-ray spectra look different. 
They show strong dependence on the observation angle in 
the GeV energy range but relatively weak dependence in 
TeV energy range (e.g. Fig. |HJ>). 7-ray spectra produced at 
the apastron passage, but at distances from the base of the 
jet of the order of a few stellar radii, do not depend strongly 
on z (see Figs.^ and d), since electrons injected at such dis- 
tance still propagate in quite similar radiation field. In fact, 
for the apastron passage of the compact object in LSI +61° 
303 (at 9.15r*), the distance of the injection place along the 
jet (z = 3r*) is located only at the distance of 9.63r* from 
the massive star (i.e very similar to the apastron distance). 
Therefore, TeV 7-ray fluxes observed from LSI +61° 303 at 
the apastron passage should be relatively larger in respect 
to their GeV fluxes in comparison to the expectations from 
the massive binary LS 5039 in which the TeV fluxes should 
be significantly lower than their corresponding GeV fluxes. 



5 PHASE DEPENDENT GAMMA-RAY 
SPECTRA AND LIGHT CURVES 

The observer is located differently in respect to the orbital 
plane of the binary systems LS 5039 and LSI +61° 303. The 
inclination angles of the binary systems are probably quite 
similar, i.e 25° in the case of LS 5039 (assuming that the 



compact object is a black hole) and 30° applied in the case of 
LSI +61° 303. However the azimuthal angles of the observer, 
measured in respect to the periastron passage, are 225° in LS 
5039 and 70° in LSI +61° 303. Therefore, phase dependent 
7-ray spectra and 7-ray light curves expected from these 
two binary systems should have different features. Below 
we show the 7-ray light curves from these binaries at the 
GeV energies (1 — 10 GeV) and the TeV energies (> 100 
GeV), based on the calculations of the IC e ± pair cascades 
occurring inside the radiation field of their massive stars. 
The phase dependent spectral features of the 7-ray emission 
are also discussed. 



5.1 LS 5039 

We have calculated the 7-ray luminosities escaping toward 
the observer as a function of the phase of the injection place 
of primary particles in the jet launched from the compact 
object in LS 5039. The light curves are obtained for two 
energy ranges, the GeV range (1 — 10 GeV) and the TeV 
range (> 100 GeV), in the case of isotropic injection of pri- 
mary 7-rays (Fig. and c) and primary electrons (Fig. [U> 
and d), with the power law spectra and spectral index 2 (as 
discussed in Sect. 4.1). We consider the injection regions at 
the base of the jet, z = (full curves) and at the distance 
2 = 5r* (dashed curves) along the jet and two inclination 
angles of the observer 9 = 25° (Fig. and b) and 6 = 60° 
(Fig.|5J; and d). In order not to complicate the geometry too 
much, only the case of perpendicular propagation of the jet 
in respect to the plane of the binary system is considered. 
The calculations of the 7-ray spectra for the jets aligned at 
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Figure 9. The power in 7-rays which escape from the binary system LS 5039 toward the observer (the 7-ray light curve) at the inclination 
angle 9 = 25° ((a) and (b)) and 6 = 60° ((c) and (d)), for the case of isotropic injection of primary 7-rays ((a) and (c)) and primary 
electrons ((b) and (d)), as a function of the phase of the injection place of the primary particles (phase is the time measured from the 
periastron passage divided by the orbital period of the binary system) at photon energies in the range 1-10 GeV (thin curves) and above 
100 GeV (thick curves). The injection place of primary particles occurs in the jet. The light curves are shown for the injection place at 
the base of the jet z m Or* (full curves) and at the distance 5r* along the jet (dashed curve). 



some angle to the plane of the binary system are straight- 
forward. They will be discussed in another work. 

If the primary particles are injected at the base of the 
jet, then the 7-ray power in the GeV energy range vary by 
a factor of ~ 2 and in TeV energy range by a factor greater 
than ~ 10 for 9 — 25° (see Fig.|5^ and b). The maximum in 
TeV 7-ray light curve occurs at the phase ~ 0.6 and corre- 
sponds to the minimum in the GeV 7-ray light curve. There- 
fore, based on these cascade calculations we predict clear 
anticorrelation between the GeV and TeV 7-ray fluxes. The 
level of variability of the GeV fluxes is larger in the case of 
injection of primary electrons in comparison to the injection 
of primary 7-rays since the primary (unabsorbed) 7-rays also 
contribute to the predicted total 7-ray light curve. There- 
fore, GeV 7-ray light curve is quite smooth for the case of 
injection of primary 7-rays. This feature might give some 
insight into the production mechanism of primary particles 
(7-rays or electrons) allowing to distinguish which particles 
are accelerated, hadrons (responsible for the primary 7-rays) 
or electrons?. If primary particles are injected farther from 
the base of the jet, then the amplitude of the 7-ray emission 
drops in both energy ranges (see dashed curves in Fig. |5Ji 
and b). However, the level of variability is still larger in the 
case of injection of primary electrons. 

The 7-ray light curves for the observer located at larger 
inclination angle, i.e. 9 = 60°, show significantly larger vari- 
ability than for 6 = 25°. For example, the 7-ray power can 
change by almost an order of magnitude in the GeV en- 
ergy range and by two orders of magnitudes in the TeV en- 



ergy range for the case of injection of primary electrons (see 
Fig.|nji). Therefore, in principle the level of 7-ray variability 
might help to distinguish between different interpretations 
of the observational data concerning the inclination of the 
binary system in LS 5039 and put some insight on the nature 
of the compact object. 

In Fig. 1101 and 1111 we show the 7-ray spectra which 
escape to the observer for a few selected phases (time from 
the periastron divided by the orbital period) of the compact 
object in all considered above scenarios for the case of the 
inclination of the binary system equal to 9 — 25° and 60°, 
respectively. The cascade 7-ray spectra produced by primary 
particles injected at the base of the jet drop suddenly above 
a few tens of GeV, reach the minimum at a few 100 GeV and 
becomes flat (spectral index <2) at higher energies. Between 
~ 0.1 — 1 TeV, the spectral index does not change signifi- 
cantly with the phase of the injection place. The TeV 7-ray 
deficit is larger for the injection of primary electrons and for 
larger inclination angles. However, the spectral indexes in 
the GeV and TeV energies are close to two in spite of such 
large difference in the level of emission. These features are 
generally consistent with the observations of LS 5039 in the 
GeV and TeV energy ranges. The spectra look completely 
different for the injection place farther along the jet. For ex- 
ample, already at the distance of 5r* from the base of the 
jet and 9 = 25°, the total cascade 7-ray spectra are almost 
independent on the phase of the compact object in the case 
of injection of primary 7-rays (see Fig. 110b ) . However, for 
larger inclination angles of the binary, e.g 8 = 60°, small ab- 
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Figure 10. 7-ray spectra which escape from the binary system LS 5039 toward the observer at the inclination angle 9 = 25° for different 
locations of the injection place of the source of primary particles ((a) and (b) - 7-rays, and (c) and (d) - electrons) defined by the 
azimuthal angle measured from the periastron: w = 0° (full curve), 90° (dashed), 180° (dot-dashed), and 270° (dotted), and at the 
distance from the base of the jet z Or* ((a) and (c)) and z = 5r* (figures (b) and (d)). 
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sorption feature starts to appear at TeV energies when the 
compact object is at the periastron. In the case of injection 
of primary electrons, the TeV 7-ray spectra are on similar 
level but differ significantly in the GeV energies. This GeV 
spectrum is flatter than in the case of injection of primary 
7-rays (spectral index close to ~ 1.5 versus a slightly flatter 
than 2), due to the lack of domination of the GeV spectrum 
by the primary 7-rays. These weak dependences of the 7-ray 
spectra in the TeV energies are due to small differences in 
the angles between the observer and the injection place of 
the primary particles farther from the base of the jet (e.g. 
considered here z = 5r*) for different phases of the binary 
system. Note however, that these angles show stronger dif- 
ferences for larger inclination angles of the binary system 
(compare Fig. HUIand Fig. lllH . 



5.2 LSI +61° 303 

The features of the 7-ray light curves and phase dependent 
spectra expected from the binary system LSI +61° 303 are 
qualitatively similar to these ones discussed above for the bi- 
nary system LS 5039 (see Figs. I12| and I13H . However, there 
also some significant differences. For example, the maximum 
in TeV 7-ray light curve of LSI +61° 303 is broader and ap- 
pears at the phase ~0.2-0.5. The level of variability at TeV 
energies is smaller (variability by a factor less than ~ 10). 
The anticorrelation between the GeV and the TeV fluxes 
is not so excellent as in LS 5039 (possible shift in phase 
by ~0.2). Significantly smaller variability is expected in the 
GeV and TeV energy ranges, if the primary 7-rays are in- 
jected already at the distance of 3r* from the base of the 
jet. The level of variability in the GeV range is larger in re- 
spect to the TeV energy range when compared to the results 
obtained fro LS 5039 in the case of the primary electrons in- 
jected at the distance of 3r* from the base of the jet. The 
absorption dips in the TeV 7-ray energy range are shallower 
than expected in LS 5039. 

In summary, propagation effects influence in similar way 
the 7-ray spectra escaping from the binary systems LSI +61° 
303 and LS 5039 when the injection place of the primary 
particles is relatively close to the massive stars (not far from 
the periastron passage and for small distances from the base 
of the jet). However, significant differences are expected for 
the injection of primary electrons farther along the jet and 
at the apastron passage as shown in Figs. I12| and ll3l 



6 CONCLUSIONS 

We have performed Monte Carlo simulations of the propa- 
gation of high energy 7-rays inside compact massive binaries 
of the microquasar type, applying as an example parameters 
of LS 5039 (recently observed in TeV 7-rays and suggested 
as a counterpart of the EGRET source) and LSI +61° 303 
(similar parameters to LS 5039 and also suggested as a coun- 
terpart of the EGRET source). In both sources the optical 
depths for TeV 7-rays are greater than unity for most of the 
propagation directions, provided that the injection place of 
primary particles (7-rays or electrons) is relatively close to 
the base of the jet launched from the inner part of the accre- 
tion disks around compact objects. Therefore, these primary 
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Figure 12. As in Fig.|5]but for the binary system LSI +61° 303 
(fi=or the inclination 9 = 30°). Primary 7-rays (figure a) and 
primary electrons (b) are injected at the distance from the base 
of the jet z ~ Or* (full curves) and z 3rv (dashed curves). 

particles initiate IC e ± pair cascades in the anisotropic ra- 
diation of the massive stars. We have calculated the 7-ray 
spectra produced in such cascades and investigated their 
basic features. It is concluded that the absorption dips in 
TeV 7-ray spectra emerging toward the observer are not so 
drastic as predicted in the earlier papers, e.g. by Bottcher 
& Dermer (2005) and Dubus (2005). This is clearly seen 
by comparing 7-ray spectra calculated with pure absorption 
(e.g Figs. and EK) with the corresponding to- 

tal 7-ray spectra produced in the cascade processes inside 
the binary system (Figs. [^Ji, [1J:, |SJ;, and[7fi). This less pro- 
nounced dips are due to the re-distribution of energy in the 
primary spectrum of 7-rays from the region above ~ 1 TeV 
to 0.1 — 1 TeV energy range. Since the soft radiation field 
created by the massive stars in these two binaries are quite 
similar during periastron passages of their compact objects, 
the 7-ray spectra produced in the cascade processes do not 
differ significantly in the case of LS 5039 and LSI +61° 303. 
However, the propagation effects are responsible for essen- 
tial differences at the apastron passage, since LSI +61° 303 
is more extended (the apastron distance in LSI +61° 303 is 
a factor of ~2 larger than in LS 5039). 

Let us concentrate at first on the features of 7-ray emis- 
sion from LS 5039. The 7-ray luminosity from this source 
observed in the GeV energy range (assuming that the iden- 
tification with the EGRET source by Paredes et al. 2000 is 
correct) is about two orders of magnitude larger than ob- 
served in TeV energy range (Aharonian et al. 2005b). How- 
ever, spectra reported in these two energy ranges are quite 
similar (spectral index close to 2). These spectral features 
can be naturally explained by the propagation effects consid- 
ered in this paper, i.e. IC e pair cascading processes, inside 
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the binary system. The primary 7-rays or electrons could be 
injected with a simple power law spectrum (and spectral in- 
dex close to 2), relatively close to the base of the jet. The 
7-ray luminosity expected in this propagation model at GeV 
energies (1 — 10 GeV) should vary by a factor of ~ 2 — 3, and 
at TeV energies (MOO GeV) by a factor of ~ 10 with the 
orbital period of the binary system LS 5039. In fact, possi- 
ble variation of the TeV signal by a factor of ~ 3 has been 
recently suggested by Casares et al. (2005), although not 
statistically significant (Aharonian et al. 2005b). Moreover, 
we show that, due to the propagation (IC e pair cascade) 
effects, the maximum in TeV 7-ray light curve should occur 
at phase ~ 0.6 (measured as an azimuthal angle from the pe- 
riastron passage, see Fig. 0. This seems to be inconsistent 
with the recent suggestion by Casares et al. (2005b) who, 
reanalyzed the HESS data with the new orbital parameters 
and, concluded that maximum of TeV emission occurs at the 
phase ~ 0.9. This discrepancy does not allow us to put any 
definitive conclusions since the error bars in TeV 7-ray light 
curve, shown by Casares et al. (2005b), are very large. How- 
ever, if real, it can give some information on the efficiency 
of acceleration process of the primary particles occurring in 
the jet (possibly linked with the efficiency of the accretion 
process) with the phase of the compact object on its orbit 
around the massive star. 

Based on the propagation calculations it is concluded 
that 7-ray spectra produced in such IC e pair cascades 
should show clear anticorrelation between the fluxes in the 
GeV and TeV energy ranges (see the light curves in Fig- - 
We also investigated the dependence of the propagation ef- 
fects on the distance of the injection place from the base of 
the jet. If primary particles are injected already at the dis- 
tance of a few stellar radii from the base of the jet (the case 



of z = 5r* is discussed), then 7-ray fluxes expected from 
LS 5039 do not vary strong enough to explain relative lumi- 
nosities in the GeV and TeV energy ranges. Either a more 
complex shape for the spectrum of primary particles injected 
into the jet is required, e.g composed of at least three dif- 
ferent power laws, or two different populations of primary 
particles are needed, or the origin of different parts of 7-ray 
spectra in different regions of the jet has to be postulated. 
Moreover, the shape of the 7-ray spectra escaping toward 
the observer depends also on the inclination angle of the bi- 
nary system which is not at present well known (estimated 
on 25° for the case of a solar mass black hole and on 60° 
for the case of a neutron star). The level of variability of the 
TeV emission for the injection place of primary particles at 
the base of the jet is clearly larger for the inclination angle 
9 = 60° than for 9 = 25°. However these differences seems to 
be to low in order to put constraints on the inclination an- 
gle of binary system LS 5039 based only on the observations 
with the present sensitivity Cherenkov telescopes. 

If the identification of the binary system LSI +61° 303 
with the EGRET source 3EG J0241+6103 is correct (see 
introduction), then average 7-ray luminosity (> 100 MeV), 
~ 8 x 10 34 erg s _1 , is only about a factor of three lower than 
in the case of LS 5039. This suggests that 7-ray properties of 
these two sources should be quite similar. Our simulations of 
the propagation of primary 7-rays show that absorption of 
7-rays is similar in both binaries at the periastron passage of 
the compact object (if primary particles are injected at the 
base of the jet) but is less important at the apastron passage 
of LSI +61° 303. The TeV fluxes from LSI +61° 303 pre- 
dicted by the propagation effects should be relatively larger 
in respect to its GeV fluxes than observed in the binary 
system LS 5039. This feature might help in detection of LSI 
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+61° 303 at TeV energies. The maximum in TeV 7-ray light 
curve predicted by the propagation effects in LSI +61° 303 
should occur at the phase ~ 0.2 — 0.4, measured from the pe- 
riastron. The TeV 7-ray signal is expected to be modulated 
by an order of magnitude with the maximum for the case of 
injection of primary electrons injected farther along the jet 
(at a distance of a few stellar radii, see the case of z = 3r» in 
Fig. I12H . Also strong variability of the GeV flux is possible 
(even by an order of magnitude) for the case of injection of 
primary electrons farther from the base of the jet (see e.g. 
the case of z = 3r* in Fig. ll2b 1. Injection of electrons at the 
base of the jet predict modulation by a factor of ~ 2 with the 
minimum when the compact object is in front of the mas- 
sive star. Re-analysis of the EGRET data indicates probable 
modulation of the GeV signal from the source toward LSI 
+61° 303 (Tavani et al. 1998, Wallace et al. 2000) with the 
maximum emission near periastron (Massi 2004). This is in- 
consistent with the predictions of studied here propagation 
effects for the case of injection of primary electrons which 
show deficit of GeV emission close to periastron. Therefore, 
if this modulation is real then the mechanism responsible 
for such modulation of GeV emission should be even much 
more efficient since considered here propagation effects work 
in opposite direction. For the case of isotropic injection of 
primary 7-rays, the GeV signal is only weakly modulated 
with the period of the binary system due to the domination 
of primary 7-rays at such low energies. 

In this paper we have only analyzed the effects of prop- 
agation of 7-rays inside these two binary systems assuming 
that efficiency of particle acceleration and spectra of injected 
particles do not depend on the location of the injection place 
of primary particles along the jet, i.e. they are uniformly dis- 
tributed along the jet, i.e. from z = Or* to at z = 5r* (LS 
5039) or 3r* (LSI +61° 303). Proper analysis of 7-ray emis- 
sion from these sources require, accept considering the IC e 
pair cascade effects, a detailed model for efficiency of accre- 
tion process, conversion of accretion energy into particles, 
acceleration of particles (their spectra), and production of 
7-rays. Such models, which unfortunately do not take into 
account the cascade effects analyzed in this paper but dis- 
cuss only production of 7-rays in the Inverse Compton pro- 
cess far away from the base of the jet, have been recently 
discussed by Bosch-Ramon, Romero & Paredes (2005) and 
Dermer & Bottcher (2005). 
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